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Abstract 

Spinning of a common aluminum automotive casting alloy A356 (Al-7Si-0.3Mg) at elevated 
temperatures has been investigated experimentally with a novel industrial-scale apparatus. 
This has permitted the implementation of a fully coupled thermomechanical finite element 
model aimed at quantifying the processing history (stress, strain, strain-rate and tempera¬ 
ture) and predicting the final geometry. The geometric predictions of this model have been 
compared directly to the geometry of the workpieces obtained experimentally. This study 
is novel in regards to both the size and shape of the component as well as the constitutive 
material representation employed. The model predictions are in reasonable agreement with 
experimental results for small deformations, but errors increase for large deformation condi¬ 
tions. The model has also enabled the characterization of the mechanical state which leads 
to a common spinning defect. Suggestions for improving the accuracy and robustness of the 
model to provide a predictive tool for industry are discussed. 
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Nomenclature 


p 

Radial penetration of roller into workpiece employed in model 

(mm) 

L 

Element edge length 

(mm) 

k c 

Thermal conductivity 

(W m -1 °C _1 ) 

c p 

Specific heat capacity 

(J kg -1 °C _1 ) 

P 

Density 

(kg m -3 ) 

a 

Coefficient of thermal expansion 

' Cc- 1 ) 

a' 

Total thermal strain 


T 0 

Reference temperature for thermal expansion 

(°C) 

Tm 

Mandrel temperature employed in model 

(°c) 

T b 

Workpiece temperature employed in model 

(°c) 


Axial/radial roller nose positions employed in model 

(mm) 

tp 

Simulated process time 

( s ) 

f3 

Taylor-Quinney factor 


Qc, f 

Heat fluxes applied to model workpiece cooling 

(W m" 2 ) 

& 1 

Maximum principal stress 

(MPa) 

°"VM 

von Mises equivalent stress 

(MPa) 

F c 

Clamping force employed in model 

(kN) 
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1. Introduction 


Often manufacturing technologies can be combined to leverage specific advantages inher¬ 
ent to each. Axisymmetric and complex shaped components can be formed through near-net 
shape casting processes alone. However, designers must account for the microstructural inho¬ 
mogeneities and potential defects that are inherent to cast components. Standard, monolithic 
forging of castings, such as through drawing or stamping, will ameliorate mechanical property 
variations, however these process have prohibitive tooling and production costs, which act 
as significant barriers to their use. As such, adopting incremental rotary forming practices, 
namely flow forming, shear forming and / or spinning, are ideal complimentary processing 
techniques for castings with axisymmetric geometries. Incremental forming is generally a 
more economical manufacturing methodology over forming process that are applied to entire 
components, such as standard forging, because forming loads and friction are minimized and 
simpler, lower cost tooling is needed. 

Even in fight of these advantages, industry has been slow to adopt incremental rotary 
forming processes owed to often costly trial-and-error approach to process development. Par¬ 
ticularly when forming at elevated temperatures, there are a myriad of process parameters: 
forming speeds, tooling geometry, material properties and history all play a role in decid¬ 
ing the final component geometry and in-service performance. Robust process models for 
industry are required to enable virtual fabrication simulation, and thereby decrease process 
development time, if the benefits of incremental forming are to be realized with any impact. 

Rotary forming processes consist of spinning an axisymmetric workpiece through a set of 
impinging tools to locally plasticize the material. Usually these tools to generate a plastic 
zone between a rigid mandrel supporting the inner workpiece diameter and a roller in contact 
with the outer diameter of the workpiece. This zone moves with the roller or tool, with 
typical applications having the material repeatedly encountering the same plastic zone to 
achieve bulk deformation. The term ‘rotary forming’ has been used to describe this class 
of incremental forming processes in the literature. Listed in order of degree of deformation 
delivered to the work piece, flow forming, tube spinning, shear forming and spinning are all 
covered by this description. For additional, more nuanced process descriptions the reader 
is directed to the literature review conducted by Wong et al. (2003) focussing on all of the 
aforementioned processes and provides detailed descriptions of each of the processes beyond 
deformation levels. More recently, Music et al. (2010) has conducted an extensive review 


focussing on the mechanics of spinning processes. 

The application of rotary forming techniques to Al-Si-Mg cast aluminum alloys has been 
studied previously by Mori et al. (2009) and Cheng et al. (2010). Both authors conducted 
spinning trials at elevated temperatures. Both studies showed an amelioration in ductility 
and slight changes in the yield strength post processing and tempering. Cheng et al. docu¬ 
mented that there was a slight improvement in corrosion resistance beyond improvement of 
mechanical properties. Mori et al. found that porosity was eliminated from A357 alloy blanks 
machined from larger castings and formed at temperatures between 350 and 400° C. Mori et 
al. also reported that mechanical properties improved during spinning at lower temperatures, 
however process-induced flaws appeared in the form of surface cracks. Their formation was 
investigated using a model developed with the commercial explicit finite element analysis, 
LS-DYNA. The sole correlation made with predictions from this isothermal process model 
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was that cracking occurred in regions with high levels of plastic strain. 

To the knowledge of the authors, there has been no previous attempt to produce a fully 
coupled thermomechanical model of a spinning operation which incorporates rate-dependent 
material properties. For aluminum alloys, the effects of strain rate and temperature have 
significant effects on the response of the material, which has not been completely captured 
in previous studies. For incremental forming operations, small regions of the workpiece will 
undergo high strain rate deformation and accompanying temperature increases due to plastic 
deformation, with the surrounding material essentially unaffected. While there are multiple 
modelling efforts of this process documented, most have employed isothermal models and 
flow stress descriptions found by quasi-static testing. While Klocke and Wehrmeister (2004) 
used different FEA packages to capture the sequential thermal and mechanical effects on 
spinning of stainless steel, this approach does not capture the heat generated due to plastic 
deformation. Furthermore, Klocke and Wehrmeister did not detail the degree of coupling 
between the thermal and mechanical aspects of their model; for example, detailing at what 
frequency the thermal and mechanical solutions were updated. Most recently, [Mohebbi 


and Akbarzadeh (2010) employed ABAQUS Explicit to generate a flow forming model for 
an annealed wrought aluminum alloy (AA6063). However, this model was only validated 
by qualitative comparisons to experimental results where forming loads were compared to 
analytical expressions derived for a different material. 


Roy and Maijer (2014) have previously employed a novel rotary forming apparatus to 


perform a series of experiments to ascertain the effects of both forming and subsequent heat 
treatment on the microstructure of A356. This study focussed on determining the effects 
of forming on precipitation and secondary phases and overall implications for further heat 
treatment. The current paper will document how these experiments have formed the basis 
of a modelling effort aimed at determining the evolution of strain, strain-rate and temper¬ 
ature during rotary forming at elevated temperatures. Roy et al. (2012) have previously 
characterised the constitutive behaviour of as-cast A356. In this latter work, a number of 
constitutive expressions to describe the material at various temperatures, strain rates and 
deformation levels were developed with varying degrees of efficacy of describing the actual 
mechanical response. This material description has been leveraged to develop the forming 
model described presently. The results of the model are compared to rotary forming experi¬ 
ments which cover a range of forming conditions and can be best described as spinning-on-air. 
In some cases, the deformation imposed was found to be more than that expected by stan¬ 
dard spinning and so as a point of convenience, the term spinning and rotary forming are 
used interchangeably. As there is a relative paucity of studies quantifying the degree to 
which standard FEA-based models of spinning compare to experimental results, the geo¬ 
metric predictions of the model are compared directly to the geometry of the workpieces 
obtained experimentally. This comparison is followed by discussion regarding proposed im¬ 
provements targeting the model’s efficacy and accuracy from the standpoint of industrial 
users as a predictive tool. 


2. Experimental forming methodology 

An experimental forming apparatus was developed to perform spinning experiments at 
elevated temperatures with workpieces representing the same geometry as an automotive 
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Figure 1: A top-down view of specific apparatus components in (a) and a parametric depiction 
with forming directions and speeds in (b). 


wheel. The apparatus consists of a belt-driven lathe driven by a 30 HP motor, a bank of 
propane torches to heat the workpiece with approximately 82 kW of total heat output and 
tooling in the form of a mandrel and roller assembly. The roller, made from AISI-8620 steel, 
had a diameter of 120 mm, a nose radius of 10 mm and was held rigidly 15° off of the radial 
axis of the workpiece (Fig. [l|. 

The mandrel was constructed from AISI-4130 and designed to account for the workpiece’s 
thermal expansion and hold the workpiece rigidly during forming, matching the workpiece’s 
ID (Fig. §. The design required that the workpiece was to be heated to forming temperatures 
from ambient conditions in situ, while held in place by manually actuated, captive clamps. 
The mandrel was tapered to ease the release of the formed workpieces, and account for 
thermal expansion of the workpiece. The mandrel was bolted directly to the spindle of the 
lathe on the non-forming end, and supported by a five center on the tailstock end, such that 
the outer diameter of the mandrel had 0.5 mm circular runout. The mandrel was furthermore 
instrumented with a series of type-K thermocouples (TCs) that were read via a wireless data 
acquisition system fixed to the quill of the lathe. In conjunction with TCs peened to a single 
commissioning workpiece, these transducers were employed to capture the thermal profile of 
the forming process. 

As-cast (AC) A356 (Al-7wt%Si-0.3wt%Mg) wheels, produced via the low-pressure die 
casting process, were provided by a North American wheel manufacturer. Workpieces or 
blanks suitable for spinning were machined from these wheels by removing the spokes and 
hub from the wheel, but preserving the main annulus. The resulting workpiece measured 
approximately 140 mm axially, had a minimum internal diameter of 330 mm and an approx¬ 
imate wall thickness of 10 mm. 
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Figure 2: Depiction of the rotary tooling assembly. 


Prior to performing forming experiments, the mandrel was preheated to approximately 
150°C with the heating system, and then sprayed with a colloidal graphite coating to act as 
a lubricant. Workpieces were preheated to the same temperature in a box furnace, and the 
inner diameter was sprayed with refractory-type coating (Foseco DYCOTE 32). This coating 
was applied to assist in workpiece removal as well as to suppress heat transfer between the 
workpiece and the mandrel, diminishing workpiece heating times. The outer diameter of the 
workpiece was coated with the same colloidal graphite as the mandrel outer diameter. The 
workpiece and tooling were then allowed to cool to ambient conditions prior to commencing 
forming trials. 


2.1. Thermal profile 

In order to characterize the thermal profile in the workpiece prior to forming, 2 type- 
K thermocouple (TC) junctions were peened into the surface of a workpiece, offset on the 
circumference by 30° at different axial locations (Fig. [3a] ) , and then the workpiece fitted 
to the mandrel. This workpiece was employed only to measure the characteristic thermal 
profile of the forming process, and was not subjected to forming. The TCs mounted on the 
workpiece were monitored in conjunction with the other 10 TCs embedded near the surface 
of the mandrel (Fig. [2] and 3a) to measure the temperature variation with time in both the 
workpiece and mandrel during preheating and forming experiments. 

Heat was applied from ambient conditions with the mandrel rotating at 20 RPM until 
thermal expansion caused the clamps to loosen, after which heating was interrupted, the 
clamps re-tightened, and heating resumed. Once the workpiece reached the forming temper¬ 
ature of ~350°C, the mandrel rotation rate was increased to the forming rate (281 RPM) to 
determine the effect of increased rotational speed on temperature. The temperatures mea¬ 
sured in both the mandrel and workpiece during this exercise are presented in Fig. |3b| and 
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(a) Thermocouple locations 




(b) Mandrel temperatures (c) Workpiece temperatures 

Figure 3: Thermocouple record and axial location on both mandrel and workpiece employed 
for thermal characterization. Solid arrow in (b) and (c) indicates start of forming. 


[3cl respectively. The dashed lines indicate where the mandrel rotation and heating system 
were disengaged to tighten the clamps, and then heating resumed. The solid arrow indicates 
the point where the mandrel rotation rate was increased from 20 to 281 RPM. 

This exercise showed that there was little temperature variation in the workpiece (within 
20°C), with the TC closest to the mandrel interface reading the coldest temperature. The 
channel exhibiting both the highest mean temperature and largest heating rate on the man¬ 
drel was located immediately below the workpiece, and the temperature and heating rate 
dropped axially. There was a 1-2°C difference in temperature between the TCs offset cir¬ 
cumferentially on the mandrel, and due to this minor difference, the temperature from these 
TCs has been presented as an average value (TM3-7 in Fig. [3j) . Increasing the mandrel speed 
from the preheating rate to the forming rate decreased the heating rate in both the mandrel 
and workpiece. The mandrel temperature continued to increase, while the workpiece temper¬ 
ature remained nearly constant at the TB2 location over twice the length of time necessary 
for forming. The temperature difference between TB1 and TB2 decreased gradually. 

These observations have been used to define the initial and boundary conditions employed 
in the model, described subsequently. 

2.2. Forming methodology 

Workpieces intended for forming did not have TCs installed owing to their interference 
with forming. Surface temperature measurements were performed manually on each work- 
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piece at the TB1 and TB2 locations with a type-K TC surface probe (Omega model number 
88108) every 3 minutes during heating, as well as immediately before and after forming. 
Workpieces were preheated to the target temperature of 355±8°C using a bank of propane 
torches. Heating times were between 17-23 minutes. The variance in workpiece heating time 
is attributed to the mandrel fitment; workpieces better conforming to the mandrel surface 
took longer to heat up as heat transfer to the mandrel was improved in spite of the refractory- 
type coating. Once the workpiece was at the forming temperature, the mandrel speed was 
increased to 281 RPM and the roller was brought into contact with the workpiece at ap¬ 
proximately 30 mm per minute. The longitudinal thread-cutting feed screw on the lathe was 
then engaged to move the roller axially at a rate of 0.21 mm per revolution (Fig. ©■ Heat¬ 
ing via the propane torches was continued during forming to lessen heat loss and maintain 
the workpiece temperature. These parameters were selected based on process information 
inferred from Cheng et^al. and the capability of the lathe. 

Once forming was complete, the workpiece was removed from the mandrel and left to air 
cool, avoiding potential distortion that could result from quenching. The deformation in all 
experiments was such that the inner diameter of the deformed region on the workpiece did 
not contact the mandrel surface. During the course of the process, the workpiece temperature 
remained above 340°C for all trials. 


2.3. Experimental results 

In this paper, the results from forming two workpieces with increasing levels of deforma¬ 
tion will be presented. The cross-sectional geometry of an undeformed workpiece as well as the 
two deformed workpieces were characterized by scanning the surface of a cross-section from 
each part and performing edge detection on the images. The profiles for each condition are 
presented in Fig. [4j The formed specimen with the least deformation (Fig. |4b| ) corresponds 
to conditions normally seen in spinning operations, where plane strain effects dominate and 
little to no workpiece wall thickness changes were observed. The highly-deformed specimen 


showed an appreciable decrease in workpiece wall thickness (Fig. 4c), which is a departure 
from a standard spinning-imposed deformation. 

The microstructure of the undeformed and deformed specimens was examined using op¬ 
tical microscopy once the surfaces of the cross-sections where polished. The deformed work- 
pieces exhibited compacted dendritic structures, as demonstrated in Fig. [5| This change in 
microstructure is similar to that reported by [Mori et al.j . It appears that grains with pri¬ 


mary dendrite arms orthogonal to the forming direction have had the secondary dendrite arm 
spacing reduced as the grain is compacted. For grains with primary dendrite arms parallel 
to the forming direction, the secondary dendrite arm spacing has been increased as the grain 
conforms to forming. 

As demonstrated in Fig. [6| the workpiece with the heaviest deformation showed surface 
defects in the form of cracking on the outer diameter in the worked regions. Cracking did 
not manifest on the low deformation sample. Cracks appeared on the surface of the sample 
early in the forming pass and alternated from opening counter to the forming direction (A, 
C) to predominately opening with the forming direction (B, D-U). The cracks do not extend 
any further than 140 pm into the bulk of the material, as measured normal to the forming 
direction. With the exception of one crack (C) extending up to 320 pm in length along the 
axis of the workpiece, the cracks tended to be short in axial length (not exceeding 200 pm). 









(a) Undeformed 




Figure 4: Cross-sectional geometries of an undeformed workpiece (a), low (b) and high (c) 
deformation workpieces. A depiction of the initial location and orientation of the roller nose 
is also provided on the cross-sections of the deformed specimens. 


The defects observed were larger in scale than those reported by|Mori et ah, who reported a 
maximum size of 60 pm and these cracks were found to predominantly penetrate the surface 
of the sample with an angle of approximately 45° from the radial direction, occurring in 
eutectic-rich regions. 

This type of localized failure is characteristic to the forming process, as it has been 
characterized in other forward rotary forming operations. Cracking or ‘fish scaling’ in rotary 
forming is caused by highly localized shear occurring both ahead of and behind the roller 
interface. The macroscopic morphology of these types of defects are described by quite well 
by Rajan and Narasimhan (2001) for ferrous material. The extent of cracking is dependent 
on processing parameters and the local strength of the material. Both of these factors decide 
the overall crack morphology. Thus, the brittle eutectic phase appears to be a weak point 
in A356 where cracks initiate. This will be discussed subsequently in conjunction with the 
results of the model. 

To eliminate cracking, either the forming temperature may be increased, or the forming 
parameters altered to avoid high levels of radial shear. Increasing the forming temperature 
from 350 to 400°C was confirmed to arrest the formation of cracks by Mori et al. for sim¬ 
ilar forming parameters. These results also agree with the study of A356 under torsion by 


McQueen et al. (1998), who reported significant increases in strain to fracture moving from 
300°C to 400°C, particularly at lower strain rates for this material. Thus, in addition to 
temperature, forming speed may be one of the process parameters that could be changed to 
reduce the frequency and severity of surface defects for rate sensitive materials. 
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(b) As-deformed 


Figure 5: Initial as-cast microstructure in (a), and resulting microstructure in (b) after 
forming. Microstructure in (a) corresponds to position (i) in Fig. 4a, and microstructure in 


(b) corresponds to position (ii) in Fig. 4c 
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(b) Micrographs 



) 200 pm 


Figure 6: Surface fracture details of heavily-deformed material. Micrographs are are aligned 
with the forming direction moving right to left. 
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3. Coupled thermomechanical model development 

In order to calculate the evolution of stress, strain, strain-rate and temperature within 
the workpiece during rotary forming, a fully coupled thermomechanical model was developed 
using the commercial FEA software package, ABAQUS. The overall model included steps 
consistent with the experimental procedure outlined above by using submodels for workpiece 
preheating, forming at elevated temperature and cooling once forming was complete. 

For the preheating model, axisymmetric heat transfer conditions were assumed to prevail 
in both the mandrel and the workpiece. Thermal and mechanical boundary conditions were 
imposed in order to calculate the preheated geometry and temperature distribution prior to 
forming. The results from this model were then translated to a 3D mesh of the workpiece. 
Forming operations applied to the workpiece were then modelled employing an adiabatic 
assumption using rigid analytical tooling to describe the mandrel and roller. The use of rigid 
analytical surfaces to describe the tooling was employed to reduce contact complications 
inherent with employing meshed instances of the tooling. Unfortunately, this assumption 
eliminates the ability to calculate the temperature history of the tooling. Mechanically, the 
choice of rigid surfaces is reasonable given the differences in flow stress and Young’s modulus 
between A356 and tool steel at the forming temperatures. The results from the forming 
model were then used to populate a final submodel where the workpiece cooling in air was 
simulated to permit direct comparison between the experimentally achieved and simulated 
geometries at ambient temperatures. 

Both implicit and explicit FEA methods were employed in simulating the overall process. 
Implicit analysis is best applied to problems without time-dependent geometric discontinu¬ 
ities as longer time steps can be used to converge on a particular solution. Meeting this 
criteria, the preheating model component used an implicit technique. However, for problems 
with large scale plasticity and involving contact, implicit methods are often unsuitable due to 
lack of convergence. In this latter case, explicit approaches are better suited as the solution is 
explicitly marched forward in time. As a result, the forming and cooling model components 
used explicit techniques. 

3.1. Geometry 

The three main components considered in the model were the workpiece, mandrel and 
roller. During preheating, only the mandrel and workpiece were considered as 2D axisym¬ 
metric descriptions. Deformation conditions present during forming necessitated a complete 
3D representation the stages occurring after preheating. 

The rigid surface geometry used for the roller was sized according to the experimental 
description in Section [2] which was assumed to retain dimensions consistent with ambient 
temperature. This assumption was made because the roller was in contact only over a small 
area for a relatively short period of time, and any change in geometry due to the small 
temperature increase would be negligible. The geometry of the rigid surface used for the 
mandrel during forming was based on the temperature-corrected geometry of the mandrel 
calculated in the preheating submodel. 

For the preheating submodel, a 2D section of the workpiece was decomposed into quadri¬ 
lateral elements with a minimum element edge length of 2 mm. The resulting mesh had 
573 nodes and 487 2D-axisymmetric, quadrilateral elements (Fig. [7J). The elements selected 
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(b) Resulting deformed mesh 

Figure 7: Preheating simulation boundary conditions and resulting deformed workpiece mesh 
used as initial conditions for forming. 


for the axisymmetric model featured hourglass control and reduced integration (CAX4RT). 
Five nodes on this 2D section were designated as tracking nodes to assist in later assigning 
boundary conditions for forming. A cross-section of the mandrel tooling was meshed with 
slightly coarser, elements having a minimum element length close to 4 mm, resulting in 673 
nodes and 515 elements. 

Four tracking nodes were designated on the surface of the mandrel to track thermal 
expansion, in addition to the five on the outer surface of the workpiece. The simulated 
displacements of these nodes during preheating were used to define the rigid analytical surface 
of the mandrel and to define where the roller contacted the workpiece axially and radially. 

To generate the mesh for use the in the forming submodel, the preheated 2D-axisymmetric 
mesh of the workpiece was revolved about its axis of symmetry to form a 3D mesh. In order 
to minimize the element count, the full 360° section of the workpiece was discretized into 
297 circumferential segments, corresponding to the maximum element edge length of the 2D 
mesh (L max ~ 4 mm, Fig. [7J). This resulted in 170181 nodes and 144639 brick/hexagonal 
elements (C3D8RT), which inherited the reduced integration and hourglass controls from 
the 2D-axisymmetric elements. This mesh density was selected primarily on the basis of 
minimizing the computational resources while ensuring at least 4 elements were in contact 
with the roller surface during forming. 

3.2. Material properties 

The model required the elastic constitutive and thermophysical properties of both the 
mandrel and workpiece, as well as a description of the rate-dependent plastic behaviour of 
A356. Material properties were assumed to be isotropic. 
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Table 1: Thermal properties of A356 and from AISI-4320. Properties for A356 are those 
employed by Hetu et al. (1998) and AISI-4320 properties are from ASM International (1990). 



C p ( J kg' 1 °C 1 ) 

k c (W m' 1 °C“ 1 ) 

a (°C 1 ) 

A356 

898.7 + 0.4270T 

7146 + 4.150T 

2.260 x 10" 7 + (2.39 x 10" 8 ) T 

AISI-4320 

452.1 + 0.4740T 

43.33- (1.090 x 10" 2 )r 

6.50 x IQ" 7 


An extended Ludwik-Hollomon constitutive expression developed by Roy et al. (2012) 
based on compression testing was implemented in the model via ABAQUS user subroutines. 
These subroutines calculate the local flow stress based on the temperature, strain and strain 
rate at each respective integration point for both the preheating (implicit) and forming (ex¬ 
plicit) simulations. This material model was developed to describe equivalent stress based on 
experimental measurements at temperatures between 30 to 500°C, strain rates between 0.1- 
10/s and plastic strain up to 0.35. The temperature-corrected elastic modulus and Poisson’s 
ratio were incorporated as tabular data based on Frost and Ashby (1982). 

The thermal-physical properties of A356 were taken from Hetu et al. (1998). Thermal 
conductivity and specific heat capacity (Ay and C p , respectively) as a function of temperature 
were implemented in the model using tabular data. The functions used to generate the data 
are shown in Table [lj A constant density, p, of 2670 kg m -3 was assumed as the density of 
A356 does not change significantly over the range of processing temperatures. The expression 
for thermal expansion, a , given by |Het u et al . and provided in Table[lJ was used to calculate 
the tangent thermal expansion data table as per the ABAQUS Analysis User’s Manual (2013). 

In order to describe the thermal expansion of the mandrel solely during the preheating 
phase of the process, thermal-physical properties for AISI-4320 were necessary. The values 
for Ay, C p , and a were entered in tabular form according to ASM International (1990) and 
are summarized in Table [lj Density was taken to be 7850 kg m -3 for AISI-4320, a modulus 
of 205 GPa and Poisson’s ratio of 0.29 were prescribed. 


3.3. Boundary and initial conditions 

Each of the simulations (preheat, forming, and cool-down) require mechanical and thermal 
boundary conditions to represent the sequence of steps taken experimentally, as previously 
described. Thermal boundary conditions were necessary to heat the workpiece and mandrel 
during preheating. During forming, the thermal boundary conditions on the workpiece were 
transformed to adiabatic conditions. Boundary conditions were then applied again to reflect 
cooling of the workpiece post-forming. Mechanical boundary conditions were used throughout 
these stages to control the relative movement of the workpiece and tooling and to describe 
the contact conditions. 


3.3.1. Preheating analysis 

In the preheating analysis, mechanical boundary conditions were applied in the form of 
node-based constraints on the spindle end of the mandrel section, such that axial movement 
was suppressed. This caused the mandrel to expand away from this constraint as the man¬ 
drel temperature increased. This constraint is consistent with the conditions present in the 
forming apparatus where the mandrel is rigidly fixed to end of the spindle and expands when 
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heated towards the tailstock via the spring-loaded centre. A mechanical contact boundary 
condition was applied between the outer diameter of the mandrel and the inner diameter of 
the workpiece. Most rotary forming simulation studies have assumed frictionless contact, e.g. 


Mohebbi and Akbarzadeh (2010) and Xue et al. (1997). A recent study conducted by Wang 


and Long (2011) on the unlubricated spinning of steel employed a Coulomb friction coefficient 


of 0.2 for tooling interactions. A Coulomb friction coefficient of 0.1 has been applied in the 
present work to all surface interactions as some friction is expected in spite of the graphite 
lubricant employed during experiments. A 10 kN force (F c in Fig. 7a) was evenly applied 


to the clamp face of the workpiece, causing the workpiece to maintain mechanical contact 
with the mandrel throughout the heating cycle and resulting thermal expansion. The load 
selected did not induce plastic deformation in the workpiece, but was sufficient to maintain 
contact. 

Attempts were made to employ inverse thermal analysis to quantify the heat flux gener¬ 
ated on the workpiece surface by the propane torches based on the experimental TC record 
(Fig. § and the domain described in Fig. [7| ft was found that the problem was severely 
ill-posed owing to the cropped mandrel domain, and the transient nature of the heat transfer 
between the workpiece and the mandrel. The variability observed in the heating rates be¬ 
tween experiments further compounded difficulties with this analysis route. The heating rate 
variability was attributed to geometric variance in the workpieces, resulting in non-uniform 
axial and circumferential contact with the mandrel. As such, uniform heating conditions in 
the form of prescribed temperature boundary conditions were applied to the nodes corre¬ 
sponding to the mandrel and workpiece in the preheating model. These boundary conditions 
were prescribed to reach the measured temperatures of T m = 220°C and T& = 350°C for the 
mandrel and workpiece, respectively, linearly ramping from an initial condition of 30°C over 


20 minutes. These boundary conditions ( dT m /dt and dT^/dt) are shown in Fig. 7a 


3.3.2. Forming analysis 

The thermal solution and geometry from the 2D preheating model was then translated 
to 3D (see Section 3.1) with the thermal solution employed as initial conditions for the 3D 
forming model. Mechanical boundary conditions for the workpiece were modified such that 


nodes on the clamp face were constrained, replacing the 10 kN force, F c (Fig. |7a| versus Fig. 

8aJ. 

Although the meshed instance of the mandrel was replaced by an analytical surface in the 
transition between preheating and forming, the mechanical contact boundary condition be¬ 
tween the workpiece and the mandrel remained unchanged. An additional contact boundary 
condition was specified between the roller and the outer diameter of the workpiece. Adhering 


to the guidelines set out by Wong et al. (2008), the workpiece was kept stationary and the 


roller was moved about and along its axis. 

During the forming experiments, the roller was moved into contact with the workpiece 
radially, and then set to move across the workpiece. In order to reflect this, the roller path 
in the forming simulation, starting with a small radial clearance from the outer diameter 
of the workpiece, initially moves radially towards the workpiece to make contact and then 
transitioned to a combined rotational and axial motion. The roller path was calculated 
based on the geometry and location of the preheated workpiece (on the mandrel) and then 
implemented in the model via a data table prescribing the motion of the reference node for 
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Figure 8: Forming model domain and boundary conditions: initial 2D axisymmetric model 
which was then translated to 3D for forming and post-forming simulations. 
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the rigid analytical surface representing the roller. 

In order to locate the roller path relative to the workpiece, the tracking nodes corre¬ 
sponding to the surface of the workpiece were used to define the geometry to match the two 
experimental forming profiles. The initial radial position of the roller nose was determined 
by linearly interpolating between the coordinates of the tracking nodes and an initial clear¬ 
ance of 0.1 mm. The initial axial position for each experimental condition was corrected for 
dilatation, and located on the geometry obtained from the preheating simulation. The path 
of the roller was then set to move first radially from the initial position u a to point u n with a 


penetration of P = 0.1 mm, while rotating about the 2 -axis (refer to Fig. 8a). Upon reaching 
u n , the roller was then moved axially to point Uf while rotating about the 2 -axis. The final 
axial position of the roller, u /, was defined as the length of the corresponding experimental 
workpiece, corrected for temperature, plus 5 mm. This 5 mm clearance was imposed to en¬ 
sure that the roller moved past the final point of contact with the workpiece at the end of the 
simulation. The proscribed axial movement of 0.21 mm per revolution and circumferential 
movement of 281 RPM of the roller matched those employed experimentally. This resulted 
in simulated forming process times of 61.93 and 87.76 s. 

As mentioned previously, forming models were run adiabatically, to characterize the heat 
developed due to dissipation or inelastic heat generation without any heat loss from the 
workpiece. A Taylor-Quinney factor of (3 = 0.9 was used. 

3.3.3. Cooling analysis 

Once the forming pass simulation was complete, the model was updated to reflect cooling 
conditions. Mechanical boundary conditions were imposed to prevent rigid body movement 
of the component during cooling, updating the conditions imposed during forming (Fig. |8a| 
versus Fig. 8c). The workpiece was then cooled down in the model by applying uniform 
surface heat fluxes to the surfaces of the workpiece (Fig. 8c). A flux of qj = 4.3 x 10 3 (T — 
25) W m -2 was applied to both the inner and outer surfaces of the workpiece and q c = 
6.4 x 10 3 (T — 25) W m -2 was applied to the clamped surface. With these fluxes applied, 
the workpiece cooled to ambient temperature in 7.5 seconds employing an explicit solution 
strategy. These fluxes were selected to limit the simulation time required to reach ambient 
temperature and to induce no further plasticity in the workpiece. Simulating the cooling 
step is necessary in order to compare the experimental workpiece dimensions at ambient 
temperatures with those predicted by the model. 

3-4- Computational overhead and scaling 

As described above, both implicit and explicit FEA methods were employed in simulat¬ 
ing the process, with explicit forming and cooling model components requiring the highest 
computational overhead. This is because the maximum explicit time step is a function of the 
speed of sound through the shortest element edge employed. 

A common approach to reduce computational resources for explicit analyses is to employ 


time or mass scaling. Both Mohebbi and Akbarzadeh (2010) as well as Wong et al. (2008) 


used these techniques to speed up solutions, but did not substantiate or justify the degree 
of scaling. This concept has been described in depth specifically for coupled FEA employing 


A356 by Roy and Maijer (2012). In order to determine the limits of scaling for the material 


model employed, Roy and Maijer used a select number of experimental thermomechanical 
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tests which were employed to populate the material model and to validate a model of these 
tests using an implicit FE formulation. These implicit results and experimental data were 
used as points of comparison to explicit simulations with increasing levels of time and mass 
scaling. This study showed that time scaling of up to 75 times could be employed before 
incoherent flow stress between that modelled implicitly and experimentally obtained data 
was observed; all forming simulations in the current study employed a conservative 50x time 
scaling factor. 

The shortest model computation times were arrived at by using parallel processing ex¬ 
tensively on a modern multi-node computational cluster. While operating with a single pro¬ 
cessor, computation times were overly prohibitive; demanding at least 30 hours per second of 
scaled process time. The effect of additional processors significantly improved computation 
time, to a certain extent. Moving to two processors showed a 104% improvement, but the 
rate of improvement showed diminishing returns with additional processors; moving from 7 
to 8 processors only displayed a 2% improvement with the additional processor. This is due 
to increased communication overhead as ABAQUS Explicit holds all element calculations on 
a single processor, and passes node-based calculations on to each additional processor. While 
adaptive meshing would be another method to improve solution times by maximizing ele¬ 
ment edge lengths and thereby maximizing time steps, ABAQUS does not support adaptive 
meshing on explicit coupled thermomechanical analyses. 

The net result is that ABAQUS requires purely Langrangian 3D domains for fully cou¬ 
pled thermomechanical simulations, thereby demanding significant computational resources. 
Even with high levels of parallelization and without adaptive meshing, more than 3 hours 
of computation time were required per second of scaled process time. This corresponds to 
approximately If days for longest running model described above. 

4. Results and discussion 

Following the development of the model, it was run to infer the forming conditions 
achieved in the experiments. This model generated a large amount of information regarding 
the process including the ability to track the evolution of the stress state, strain rate, and 
temperature in the workpiece during forming. The results of the forming model reflecting the 
experimentally formed workpieces will be presented as an example of the model capabilities. 
These results will be followed with a comparison of the final workpiece geometries generated 
experimentally and those predicted with the model. Finally, the stress state imposed on the 
workpiece during high-level deformation will be presented and discussed within the context 
of surface defect formation. 

4-1. Model Results 

The forming model provides a great deal of insight into the overall process. One facet 
is the evolution of the stress-state and deformation of the overall workpiece during forming. 
This is demonstrated in Fig. [9| which shows the distribution of the equivalent stress on the 
surface of the least-deformed workpiece at the start, midpoint and end of forming. This 
corresponds to three roller positions: i) u n , ii) at the midpoint of the roller travel, and iii) 
just prior to the roller leaving the surface of the workpiece. In this figure, the workpiece has 
been rotated such that axis of symmetry of the roller is parallel with the centerline of the 
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workpiece, affording an orthogonal view of the workpiece along its axis during forming. Here, 
and further, von Mises equivalent stress (ct V m) has been used to analyze the workpiece stress 
state owing to the multiaxial nature of loading. 


Overall, the distribution of surface stresses are similar to those demonstrated by Xu et al. 


(2001) found by simulating the spinning of tubes: a very high compressive region in proximity 
to the roller acting radially, with a tensile region outside of this. However, due to the irregular 
shape of the workpiece employed, the forming zone is highly dynamic as compared to tube 
spinning. 

At the start of forming, the forming zone is highly localized and close to being symmetric 
about the centerline of the workpiece. The majority of this zone, as identified by the stressed 
regions greater than 125 MPa, is directed circumferentially along the path of the roller, 
with the maximum contact stress appearing slightly behind the centerline, relative to the 
workpiece rotation. At the midpoint of forming, the stress state has evolved such that the 
majority of the forming zone remains ahead of the roller both axially and circumferentially, 
and significant stresses have evolved elsewhere on the surface. The region ahead of the roller 
and towards to the end of the workpiece along 2 shows significant load directed axially. Also 
at this stage in forming, the overall bounds of the circumferential region carrying an elevated 
load has increased dramatically from the start of forming. At the end of forming, the stress 
state returns to being highly localized immediately beneath the roller, behind the centerline 
relative to the workpiece rotation. 

Fig. [lO] shows an oblique view of the workpiece at the same stages during forming pre¬ 
sented in Fig. [9| Also shown in this figure are the nodes in contact with the roller at each 
stage. At the start of forming, a faint ridge can be seen due to the initial contact of the 
roller, with a relatively small roller contact footprint. At the midpoint of forming, this ridge 
is more visible and contact has extended to include a small pileup of material ahead of the 
roller, which results in a larger effective contact region. The region ahead of the roller in the 
axial direction along z demonstrates slight diametral growth, coinciding with the high axial 
stress levels as compared to behind the roller (along z). At the end of forming, the contact 
patch has diminished according with the dissipation of material pileup. In this last stage, the 
formed regions behind the roller exhibit ridges attributed to non-uniform pileup dissipation 
as the roller progressed across the surface of the component. 

Orthogonal views of the deformation and stress states occurring in the highly deformed 
workpiece at the 3 stages during forming are presented in Fig. [Tl} Overall, the stress 
magnitude is significantly higher than seen with the least-deformed workpiece, coinciding 
with a more aggressive forming profile. Similar to the other workpiece, the peak stresses on 
the surface of the workpiece occur slightly ahead of the centerline at all forming stages. At the 
start of forming, the stress state is similar to the least-deformed workpiece in that the contact 
stress is approximately axially symmetric about the centerline of the roller. However, this 
stress state embodies a much larger region projected towards the fixed end of the workpiece, 
which is attributed to bending stresses caused by the roller contact. Midway through forming, 
the stress state has changed dramatically to be quite disparate from the least-deformed 
workpiece. This is attributed to the workpiece buckling and forming a convex flange ahead 
of the roller, which was observed experimentally. At this stage in forming, two regions of 
elevated stress appear on the workpiece circumference. One aligned with the centerline of the 
roller, and the other appearing on the flange, behind the roller circumferentially and ahead 
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Figure 9: Predicted equivalent stress state (ctvm) on the surface of the least-deformed work- 
piece during forming. 



















































(a) Start 



(c) End 


Figure 10: Oblique views of the simulated equivalent stress state immediate to the roller on 
the surface of the least-deformed workpiece. Inset shows nodes in contact with the roller. 
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of the roller axially. At the end of forming, the stress state returns to being localized to the 
vicinity of the roller, however, the axial length shows a large degree of irregularity around 
the circumference, which was not seen in the least-deformed workpiece. 

Oblique views of the highly-deformed workpiece at the same forming stages are shown 
in Fig. [12] in conjunction with the nodes in contact. This shows that the same number 
of elements are in contact with the roller as the least-deformed workpiece at comparable 
forming stages. However, this is the sole similarity. The initial ridge formed by roller contact 
at the start of forming is much more pronounced. Beyond the formation of the flange midway 
through forming, the surface previously encountering the roller is much more irregular and 
a small region of pileup is seen at the edge of the flange closest to the roller. At the end of 
forming, the flange has collapsed, and the edge of the workpiece shows localized irregularities. 
The ridges seen in the least-deformed workpiece that were attributed to non-uniform pileup 
dissipation are much more exaggerated in this forming case. However, the ridges in the 
highly-deformed case are much less radially consistent than in the least-deformed workpiece. 
Clearly, the simulation predicts significantly less uniform deformation than observed in the 
least-deformed workpiece, with the workpiece wrinkling as the flange buckles. 

Two orthogonal views of the highly-deformed workpiece at the last stage of forming are 
provided in Fig. [13] with contours showing the predicted axial displacement, u z . This result 
highlights the non-uniform distribution of axial deformation, in particular the formation of 
lobes that are apparent at the end of forming. This phenomena is discussed further in the 
next section. To further examine the evolution of other process variables beyond stress and 
displacement, cross-sectional views of the simulated workpieces immediately underneath the 
roller have been extracted at each stage of the forming process. An example of one of these 


locations is shown in Fig. 13 


Contours of equivalent stress, equivalent plastic strain, strain rate, and temperature on 
the cross-sectional planes at each of the three stages during forming of the least-deformed 
workpiece are shown in Figs. [14] - EE3 respectively. The equivalent stress in the least- 
deformed workpiece (Fig. 0 at the start of forming shows that the region of elevated stress 
is localized directly beneath the roller and does not extend very far through-thickness. As 
forming progresses, the stress distribution evolves ahead of the roller, staying predominantly 
localized to the outer diameter of the workpiece as it bends to conform with the forming 
profile. The distribution of plastic strain (Fig. 15) mirrors the stress profile, with a peak 
of 0.7 on the outer diameter and 0.3 on the inner. The strain rate distribution (Fig. 0 
which encompasses both elastic and plastic strain rates, has a peak located at the roller 
interface at the start of forming. Midway through forming, the peak strain rate has shifted 
to being slightly ahead of the roller and midway through the workpiece thickness. This is 
attributed to the combined effects of surface deformation and bending. The peak rate remains 
at approximately 4-4.5 s -1 for most of the forming pass, rising to 7 s -1 briefly at the end. 
Reflecting the relatively low amount of strain imparted to the workpiece, the temperature 
has only increased by approximately 10°C (Fig. 0. less than 3% of the initial temperature. 
This implies that for these forming conditions a small amount of heat was generated and a 
fully coupled model was not necessary (i.e. conditions were near isothermal). 

Contours of equivalent stress, equivalent plastic strain, strain rate, and temperature on the 
cross-sectional planes at each of the three stages in forming of the heavily-deformed workpiece 
are shown in Figs. [18] - (2ll respectively. The equivalent stress state occurring radially in 
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(c) End 

Figure 11: Simulated equivalent 
workpiece during forming. 


(b) Midpoint 

stress state (ct V m) on the 


(a) Start 

surface of the highly-deformed 
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(b) Midpoint 



(c) End 

Figure 12: Oblique views of the simulated equivalent stress state immediate to the roller on 
the surface of the highly-deformed workpiece. Inset shows nodes in contact with the roller. 


24 













(a) Axial view (b) Normal to axis 


Figure 13: Contour plots of nodal displacement along the heavily-deformed workpiece axis, 


corresponding to Fig. 11c and 12c 
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(a) Start 



50 100 150 200 

(c) End 


Figure 14: Equivalent stress distribution on 
cross-sections of the least-deformed workpiece 
at different forming stages. 


(a) Start 



0 0.2 0.4 0.6 

(c) End 


Figure 15: Equivalent plastic strain distribu¬ 
tion on cross-sections of the least-deformed 
workpiece at different forming stages. 
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(a) Start 



(c) End 

Figure 16: Strain rate distribution on cross- 
sections of the least-deformed workpiece at 
different forming stages. 


(a) Start 



350 355 360 

(c) End 


Figure 17: Temperature distribution on cross- 
sections of the least-deformed workpiece at 
different forming stages. 
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50 150 250 

(c) End 


Figure 18: Equivalent stress distribution on 
cross-sections of the highly-deformed work- 
piece at different forming stages. 



(b) Midpoint 



0 1 2 

(c) End 


Figure 19: Equivalent plastic strain distribu¬ 
tion on cross-sections of the highly-deformed 
workpiece at different forming stages. 


the highly-deformed workpiece (Fig. 18) demonstrates a much higher load both on the outer 
and inner diameter than the least-deformed workpiece. This is likely due to bending stresses 
developed as soon as the roller contacts the workpiece. As forming progresses, the stress state 
evolves ahead and to a lesser extent, behind the roller. The highest stresses appear in the 
middle of the workpiece, differing from the least-deformed workpiece as the deformation in the 
highly-deformed workpiece is dominated by bending. The distribution of equivalent plastic 
strain (Fig. |I9| ) is primarily localized on the outer and inner diameter of the workpiece at all 
forming stages. This is similar to the least-deformed workpiece, albeit peak strains are an 
order of magnitude higher. Additionally, at the midpoint of forming, an appreciable plastic 
zone has developed on the edge of the flange, well ahead of the roller, which is attributed 
to the start of buckling in the flange. At the end of forming, the flange region exhibits the 
largest amount of strain caused by the combination of buckling and roller contact during 
forming. The strain rate (Fig. 20) distribution is similar to that observed in the least- 


deformed workpiece, however, the magnitude is approximately 4 times greater owing to the 
large amount of buckling. As there was significantly higher amounts of strain applied, the 
temperature increase is significantly higher, with the peak temperature increasing by 40°C 
(Fig. 0. This rise in temperature has a significant effect on the flow stress of the material, 
and therefore a coupled model is necessary to accurately capture the material behaviour. 


4-2. Geometric Comparison to Experimental Results 

In order to gauge the effectiveness of the model to predict the final shape of the workpiece, 
the results of the model after cooling the deformed workpiece to room temperature have been 
compared in two ways. First, the predicted workpiece lengths are compared over the entirety 
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(c) End 

Figure 20: Strain rate distribution on cross- 
sections of the highly-deformed workpiece at 
different forming stages. 



350 370 390 

(c) End 

Figure 21: Temperature distribution on cross- 
sections of the highly-deformed workpiece at 
different forming stages. 


of the circumference. Second, the predicted workpiece cross-sections are compared through 


thickness with corresponding experimental profiles (Fig. [4b] and 4c). 

Fig. [22] shows variation of the predicted axial length of the work piece around the cir¬ 
cumference for the least-deformed case. The model data shows that there are four discrete 
minima and maxima appearing along the circumference of the part, corresponding to faint 
lobes. The measured extents of the least-deformed workpiece have been plotted as dashed 
lines for comparison. The solid line coincides with the overall length of the corresponding 
experimental cross-section (Fig. [4]). The experimental workpiece did not exhibit lobes as the 
minimum and maximum length were offset circumferentially by 180°. However, with the ex¬ 
perimental workpiece measuring 149.25±0.55 and the simulated 148.43±0.62, the agreement 
between the model and experiment is within the length of an undeformed element. 

Fig. 23 compares the shape of the experimental cross-section (solid outline) with that 
predicted by the model (mesh) at the location of the maxima (on lobe) and minima (off 
lobe). The model predicts the development of a convex shape axially along the outer di¬ 
ameter, whereas the experimental profile is slightly concave. Both predicted profiles exhibit 
approximately the same amount of error in describing the experimental cross-section, with 
the inner radius being approximately 4.5 mm smaller at the midpoint of the deformed region 
than found experimentally. 

Fig. [24] presents the same comparison as Fig. [22] for the heavily-deformed workpiece. 


For the high deformation conditions, the circumferential variation in length predicted by 
the model is much more pronounced, demonstrating much larger lobes. The final length of 
the experimental workpiece was 162.15±0.65 mm, where as the model predicts a length of 
155.15±3.15 mm. The model does not predict the overall part length for the highly-deformed 
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Figure 22: Final axial length of the least de¬ 
formed workpiece compared to experimental 
measurements. 


Figure 23: Cross-sectional comparison of 
model and experimental results for least- 
deformed workpiece. 


conditions nearly as well as the least deformed forming profile. This is also reflected in the 
comparison of the experimental cross-section with model predictions at the on and off lobe 
positions (Fig. 25). The on-lobe position agrees with the experimental cross section for the 
majority of the length, with the exception of the very end of the workpiece. The off-lobe 
position departs significantly, predicting an inner radius approximately 9 mm smaller at the 
midpoint of the deformed region. 

The discrepancy between the model prediction and the experimental results for the 
heavily-deformed workpiece has been attributed to two potential sources of model inade¬ 
quacy. First, there are limitations in the model to predict localized material failure. The 
model currently predicts localized deformation, but does not predict material fracture. Ex¬ 
perimentally, surface cracks were observed on the outer diameter of the workpiece (Fig. |6|. 
Incorporating the effects of surface cracking in the model would affect the bending stresses 
significantly as the overall part stiffness would decrease. This would in turn modify the flange 
buckling phenomena. Based on the experimental results, it is surmised that surface cracking 
diminishes wrinkling and limits the formation of lobes. Second, the adiabatic conditions im¬ 
posed may reflect an inaccurate account of how heat was transferred through the component 
during forming. This would induce a decreased yield point ahead of the roller, leading to 
inaccuracies in capturing the overall forming zone size, also modifying the flange buckling 
phenomena. 

However, the model does show reasonable agreement with the experimental geometry in 
the case of the least deformed workpiece. It also provides an adequate prediction of the radial 
cross-section of the highly-deformed workpiece on lobes. The lack of agreement elsewhere for 
the highly-deformed workpiece can be attributed to cracking, and potentially inadequately 
imposed thermal boundary conditions. This suggests that the the basic deformation mecha¬ 
nism has been successfully modelled for some forming conditions. 
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Figure 24: Final axial length of the heavily- 
deformed workpiece compared to experimen¬ 
tal measurements. 


Figure 25: Cross-sectional comparison of 
model and experimental results for the 
heavily-deformed workpiece. 


4-3. Surface defect formation 


As discussed previously, [Mori et al. demonstrated experimentally that surface cracks 
similar to those found in the present study occurred in regions of high levels of strain. These 
levels of strain were identified with a strain rate independent model. In reviewing the results 
of torsion testing conducted by McQueen et al. (1998), fracture in A356 was observed to 
occur at strains of 1 and 1.5 at 300 and 400°C, respectively, for strain rates up to 5 s -1 . 
The torsion tests also showed that the equivalent stress at fracture was seen to increase with 
higher strain rates. 

Based on this information, the forming model can be used to explain the lack of obvious 
surface cracks appearing in the least deformed workpiece. The peak strain predicted by the 
model for this workpiece was 30% less than the fracture strain identified by |McQueen et al., 
at approximately the same temperature and strain rate conditions. In the case of the highly- 
deformed workpiece, the model predicts equivalent plastic strains of 1.5 or more during early 
forming, prior to the flange buckling. Furthermore, the predicted strain rate is significantly 
higher than the range employed by [McQ ueen et al., which would decrease the fracture strain 
to a greater extent. However, strain and strain rate alone are not sufficient to predict local 
failure. 

Fig. 26 shows a contour plot of the the principal stress (<ti) magnitude occurring on a 
cross-section of the highly-deformed workpiece. Overlaid on the contour plot is a quiver plot 
showing the orientation of oq relative to the workpiece axis, with origins located at the do¬ 
main’s integration points. This state reflects 17.5 mm of axial roller travel, or approximately 
25% through the forming profile, coinciding with a position halfway between those shown 
in Fig. |lla| and |llb[ The region in contact with the roller demonstrates a high degree of 
compressive stress, however the element immediately behind the roller shows a nearly zero 
stress state at the surface. The stress state is increasingly tensile moving towards the inner 
diameter. In Fig. [6j the crack morphology in terms of opening angle from the surface of the 
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Figure 26: Simulated <j\ magnitude and orientation from the forming direction immediate to 
the roller. The depicted state was extracted from the simulation reflecting heavy deformation, 
at a quarter way through the axial travel of the roller, on the z — R plane bisecting the roller. 


workpiece are quite similar. Notable exceptions include cracks A and C where the opening 
direction is reversed, but a much shallower crack is observed as compared to the majority. 
The most prevalent crack morphology matches the stress disparity as annotated by the black 
dashed line in Fig. [26j reflecting unconstrained plastic flow. 

Attributing the precise mechanical state shown in Fig. [26] to categorically identify the 
conditions necessary for crack development during forming is speculative at this stage in 
model development. This is because the material model employed is limited by the data 
available regarding fracture conditions. The model does, however, provide a framework to 
include further data with which forming parameters may be modified to conclusively mitigate 
this phenomena from occurring. One possible extension would be a criterion which tracks 
principal stresses to ascertain the degree of plastic constraint as a tool to predict whether or 
not cracking will occur. However, it will be necessary to experimentally develop the crack 
formation behaviour of this material at elevated temperatures. 


5. Summary 

In order to investigate spinning of full-size A356 components at elevated temperatures, an 
experimental forming apparatus was constructed. Two experimentally deformed components 
with different levels of spinning and buckling deformation were produced with this apparatus 
and were analyzed in this study. The workpiece with higher levels of deformation was found 
to exhibit surface defects in the form of cracking or ‘fish-scaling’. 

The various thermal and mechanical processing steps involved in processing were modelled 
as a coupled thermomechanical process in ABAQUS. This consisted of an implicit FE sub¬ 
model for preheating the workpiece to forming temperatures, as well as explicit FE submodels 
of the forming operation at elevated temperatures and final cooling to room temperature. 
The overall modelling effort contributes significantly to the overall understanding of rotary 
forming as a whole; heretofore, strain rate and temperature dependency in modelling efforts 
have not been considered in rotary forming, much less for A356. 

The model shows that the local effects of the roller interface dominate throughout form¬ 
ing, and that the strain rates achieved are quite high, approaching that of standard metal 
turning. Changes in material state, manifesting with the evolution of strain rate and tem- 
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perature produce non-uniform deformation in the form of wrinkling and the formation of 
lobes. Comparing the geometry of the lightly deformed workpiece with that predicted by the 
model shows reasonable agreement with the experimental geometry, demonstrating that the 
model successfully describes the basic deformation mechanism. Comparing the results for 
the highly-deformed workpiece showed that the model did not predict the geometry nearly 
as well. This may very well be due to the inability of the model to predict local damage. As 
additional data on cracking of A356 at high strain rates and strains at elevated temperatures 
becomes available, this model may provide the ability to modify forming parameters such 
that rotary forming defects may be avoided for this material. 

As in a number of other previous studies which have attempted to model incremental 
rotary processes using commercial code, prohibitive computational requirements remain the 
primary limitation. This is particularly true when considering thermomechanical coupling. In 
order for FEA process modelling of rotary forming to be adopted by industrial practitioners 
for full-sized components, specialized code and/or custom hardware may be required. This is 
essential if modern techniques such as adaptive meshing are to be employed and to provide 
accurate results in a reasonable amount of time. 
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